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Abstract. In this paper we finish the intensive study of three-dimensional 
Dirichlet stereohedra started by the second author and D. Bochis, who showed 
that they cannot have more than 80 facets, except perhaps for crystallographic 
space groups in the cubic system. 

Taking advantage of the recent, simpler classification of three-dimensional 
crystallographic groups by Conway, Delgado-Friedrichs, Huson and Thurston, 
in a previous paper we proved that Dirichlet stereohedra for any of the 27 "full" 
cubic groups cannot have more than 25 facets. Here we study the remaining 
"quarter" cubic groups. With a computer-assisted method, our main result 
is that Dirichlet stereohedra for the 8 quarter groups, hence for all three- 
dimensional crystallographic groups, cannot have more than 92 facets. 



1. Introduction 

This is the last in a series of four papers (see [TJ [5J [T7]) devoted to bounds on 
the number of facets that Dirichlet stereohedra in Euclidean 3-space can have. 

A stereohedron is any bounded convex polyhedron which tiles the space by the 
action of some crystallographic group. A Dirichlet stereohedron for a certain crys- 
tallographic group G is the Voronoi region Vor Gpip) of a point p € R 3 in the Voronoi 
diagram of an orbit Gp. 

The study of the maximum number of facets for stereohedra is related to Hilbert's 
18th problem, "Building up the space with congruent polyhedra" (see |12M16j ). Bie- 
bcrbach (1910) and Reinhardt (1932) answered completly the first two of Hilbert's 
specific questions, but other problems related to monohedral tessellations (i. e., 
tessellations whose tiles are congruent) remain open. An exhaustive account of 
this topic appeared in a survey article by Grunbaum and Shephard [TT] , where our 
problem, to determine the maximum number of facets — or, at least, a "good" upper 
bound — for Dirichlet stereohedra in R 3 , is mentioned as an important one. Previous 
results on this problem are: 

• The fundamental theorem of stereohedra (Dclaunay, 1961 [B]) asserts that 
a stereohedron of dimension d for a crystallographic group G with a aspects 
cannot have more than 2 d (a + 1) — 2 facets. The number of aspects of a 
crystallographic group G is the index of its translational subgroup. Delone's 
bound for three-dimensional groups, which have up to 48 aspects, is 390 
facets. 

• The three-dimensional stereohedron with the maximum number of facets 
known so far was found in 1980 by P. Engel (see [7] and |11[ p. 964]). It is 
a Dirichlet stereohedron with 38 facets, for the cubic group I4i32, with 24 
aspects. 
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(1) Bounds after processing triad rotations. 

(2) Bounds after diad rotations with axes parallel to the coordinate axes. 

(3) Bounds after diagonal diad rotations. 

(4) Bounds after intersecting with planar projections. 

Table 1 . Bounds for the number of facets of Dirichlet stereohedra 
of quarter cubic groups 



There is agreement among the experts (see [7J page 214], [TTJ page 960], [TH1 
page 50]) that Engel's sterohedron is much closer than Delone's upper bound to 
having the maximum possible number of facets. Our results confirm this. 

In 2000, the second author and D. Bochi§ gave upper bounds for the number of 
facets of Dirichlet stereohedra. They did this by dividing the 219 affinc conjugacy 
classes of three-dimensional crystallographic groups into three blocks, and using 
different tools for each. Their main results are: 

• Within the 100 crystallographic groups which contain reflection planes, the 
exact maximum number of facets is 18 [I]. 

• Within the 97 non-cubic crystallographic groups without reflection planes, 
they found Dirichlet stereohedra with 32 facets and proved that no one can 
have more than 80. Moreover, they got upper bounds of 50 and 38 for all 
but, respectively, 9 and 21 of the groups 

• They also considered cubic groups, but they were only able to prove an 
upper bound of 162 facets for them [3]. 

In [17] we improved the bound for 14 of the 22 cubic groups without reflections 
planes, the 14 "full groups" : 

Theorem 1.1. Dirichlet stereohedra for full cubic groups cannot have more than 
25 facets. 

In this paper we give an upper bound for the remaining cubic groups: the 8 
"quarter groups". It has to be noted that to get these bounds, contrary to the 
ones in the previous papers of this series, computers are used. The upper bound 
we obtain for each quarter group is shown in Table [TJ Columns (1) to (4) are the 
bounds obtained in different phases or our method, the column labeled "Final" is 
our final bound. Globally, we get the following. 

Theorem 1.2. Dirichlet stereohedra for quarter cubic groups cannot have more 
than 92 facets. 

For the sake of completeness, we include here the full list of other crystallo- 
graphic groups for which the bounds proved in this series of papers is bigger than 
38 (Table [2]). This list is the same as Table 2 in [2]. 
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Table 2. Non-cubic groups where our upper bound is larger than 38 



Theorem 1.3. Three dimensional Dirichlet stereohedra cannot have more than 92 
facets. They can possibly have more than 38 facets only in one of the 29 groups 
listed in tables[l\ and[^ 



2. Preliminaries and outline 

2.1. "Full" and "quarter" cubic groups. Our division of cubic groups into 
"full" and "quarter" ones comes from the recent classification of three-dimensional 
crystallographic groups developed in [5] by Conway et al. They divide crystallo- 
graphic groups into "reducible" and "irreducible" , were irreducible groups are those 
that do not have any invariant direction. It turns out that they coincide with the 
cubic groups of the classical classification. Conway ct al. define odd subgroup of an 
irreducible group G as the one generated by the rotations of order three, and show 
that: 

Theorem 2.1 (Conway et al. [5]). (1) There are only two possible odd sub- 
groups of cubic groups, that we denote F and Q. 
(2) Both F and Q are normal in Isom(R 3 ). Hence, every cubic group lies 
between its odd subgroup and the normalizer Af(F) and N{Q) of it. 

The second property reduces the enumeration of cubic space groups to the enu- 
meration, up to conjugacy, of subgroups of the two finite groups JV(F)/F and 
Af(Q)/Q. N(Q)/Q is dihedral of order 8 and N(F)/F has order 16 and contains 
a dihedral subgroup of index 2. 

The main difference between F and Q is that Q only contains triad rotations 
whose axes are mutually disjoint, while some triad rotation axes in F intersect one 
another. Q is a subgroup of F of order four and because of that Conway et al. call 
full groups those with odd subgroup equal to F and quarter groups those with odd 
subgroup equal to Q: Quarter groups contain only a quarter of the possible rotation 
axes. There are 27 full groups (14 of them without reflection planes) and 8 quarter 
groups (none with reflection planes, because Af(Q) does not contain reflections). 

2.2. The structure of quarter cubic groups. Throughout the paper we use 
the International Crystallographic Notation for three-dimensional crystallographic 
groups; see, e. g., [IB] . 

All the quarter cubic groups have the same odd group Q, a group of type P2i3. 
By definition, Q is generated by triad rotations. More precisely, through each point 
(x,y,z) £ (Z/2) 3 exactly one rotation axis passes, with vector: 
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Figure 1. The body-centered cubic lattice /. The triad rotations 
(in grey) together with translations of length two in the coordinate 
directions generate Q 



(1, 1, 1) if £ = y = z (mod 1) 

(—1, 1, 1) if y ^ x = z (mod 1) 

(l,-l,l)ifz^x = y(mod 1) 

(1, 1, —1) if x ^ y = z (mod 1) 
In other words, exactly one of the four diagonals of each primitive cubic cell of 
the lattice (Z/2) 3 is a rotation axis (see Fig[I}. The translational subgroup of Q is a 
cubic primitive lattice with vectors of length one. Af{Q) is the group of symmetries 
of the set of triad rotation axes. It is a group of type J-j3|, and its translational 

subgroup is a body centered lattice generated by the vectors (±^, ±2, ±^)- 

Af(Q)/Q has order 8. In fact, it is isomorphic to the dihedral group D$. Hence 
there are 8 groups between Af(Q) and Q, including both. Their lattice is drawn in 
Figure [5] In Figure [5J we represent them graphically with the conventions of [T3] : 
for each group, the intersection of a generic orbit with the primitive cell [0, l] 3 is 
considered. Because all the groups contain the triad rotation on the diagonal axis 
x = y = z, only a third of the orbit is necessary in order to describe the group. E. g., 
the orbit with respect to the subgroup that sends horizontal planes to horizontal 
planes. This is what is drawn in the figures, projected over the XY plane. 

The full dot near the origin represents a base point (a, b, c) with < a,b,c < 1/4. 
Next to each of the other orbit points there is a number h 6 {0,1/4,1/2,3/4}, 
omitted whenever it equals 0. If an orbit point is represented as a "full dot", then 
its z-coordinate equals c + h. If it is represented as an "empty dot" , then it equals 
h — c (or 1 — c, if h is zero). For extra understanding of the group, each dot is 
drawn dark or light depending on whether it is obtained from the base point by an 
orientation preserving or an orientation reversing isomctry in G. 

2.3. Outline of our method. Let G be one of the eight quarter groups and let 
p S M 3 be a base point for an orbit Gp, so that the Dirichlct stercohedron we want 
to study is the closed Voronoi region Vor Gp(p)- There is no loss of generality in 
assuming that p has trivial stabilizer, since Koch JT3J H31 completely classified cubic 
orbits with non-trivial stabilizer (more generally, those with less than 3 degrees of 
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Figure 2. The eight quarter groups. 

freedom) and showed that they produce Dirichlet stereohedra with at most 23 
facets. 
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If p has trivial stabilizer, then a small perturbation of its coordinates can only 
increase the number of facets [T7J Lemma 3.1], so we assume p to be sufficiently 
generic. 

Our method is as follows: We consider a certain tessellation T of the three- 
dimensional Euclidean space, which we call the auxiliary tessellation. For each 
tile T of T we call extended Voronoi region of T any region VorExt G T with the 
following property: 



The extended Voronoi region is not uniquely defined. Part of our work is to compute 
one that is as small as possible, for each tile T. 

Now, if To is the tile of T that contains the base point p, we call influence region 
°f To (or, of p) the union of all the tiles of T whose extended Voronoi regions meet 
VorExtc(To) in their interiors. We denote it by Infl G (To). Observe that, strictly 
speaking, Infl G (To) depends not only on G, but also on T and on our particular 
choice of extended Voronoi regions. It is easy to prove: 

Theorem 2.2. If To is the tile containing p, then all the elements of Gp that are 
neighbors of p in the Voronoi diagram Vor Gp (p) ^ e in InflcC^o)- 



That is, to get a bound on the number of facets of Dirichlct stereohedra for base 
points in To we only need to count, or bound, \Gp (~1 Infl G (To)|. 

The above "sketch of a method" is impossible to implement directly. It involves 
the computation of an infinite number of extended Voronoi regions and influence 
regions. The key property that allows us to convert this into a finite computation 
is the following easy lemma: 

Lemma 2.3. Let M{G) denote the normalizer of a crystallographic group G. IfT\ 
and T2 = pT\ are tiles related by a transformation p G M{G), then VorExt G (T2) = 
/9VorExt G (Ti) and Infl G (T 2 ) = pInfl G (Ti). 

Proof. It follows from the fact that Vor Gp ( 9 ) (p(q)) = /o(Vor Gg (g)) for every point q, 
if p is in the normalizer Af(G) of G. □ 

So, what we do is to use an auxiliary tessellation whose tiles lie in a small number 
of classes modulo the normalizer of G. Only this number of influence regions needs 
to be computed. 

2.4. Computation of extended Voronoi regions. The basic idea for our com- 
putation of extended Voronoi regions is that each translation and rotation present 
in G implies that a certain region of space can be excluded from VorExt G (To). 

Lemma 2.4. Let To be a convex domain and let p be a rotation in G with axis I 
and rotation angle a. Assume that £ does not intersect Tq. Let Hi and H2 be the 
two support half-planes of To with border in £. Let H[ and H' 2 be the half-planes 
obtained by rotating H\ and H2 "away from To" with angles ±a/2 and axis t. 

Then, for every p £ To, the Dirichlet region Vor Gp (p) is contained in the (perhaps 
non-convex) dihedral region bounded by H[ U H' 2 and containing To- 
Proof. In fact, something stronger is true: Vor Gp (p) is contained in the dihedral 
sector with axis in t, angle a, and centered at p. This follows from the fact that 
the two facets of this dihedron are parts of the bisectors of p and pp and of p and 
respectively. See the left part of Figure [3] □ 



Vor G9 ( g ) C VorExt G (T). 



Proof. See 0. 



□ 
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FIGURE 3. Situations in the Lemmas 12^1 and l2~5l 




Figure 4. Extended Voronoi region induced by three rotations 
and two translations 

Lemma 2.5. Let Tq be a convex domain and let v be the vector of some translation 
in G. Let Hi and H2 be the two support planes of Tq orthogonal to v. Let H[ and 
H' 2 be the planes obtained translating Hi and Hi "away from Tq " by the vectors 
±v/2. 

Then, for every p £ To, the Dirichlet region VoiGpip) is contained in the strip 
between H' x and H' 2 . 

Proof. Similar to the previous one. See the right part of Figure [3] □ 

So, our method for constructing extended Voronoi regions consists in identifying 
a certain number of rotations and translations in G and taking as VorExtG(Xo) 
the intersection of the regions allowed in Lemmas 12.41 and 12.51 Figure @] shows an 
example of what we mean, in which we consider three rotations and two translations 
in the plane. 

The problem with this approach is that the regions obtained in Lemma 12.41 may 
not be convex. Observe that we are going to use more than 250 rotations (see 
Tables 1^1 ITTJl [TT1 and [T2"|) . So, in principle, to compute an extended Voronoi region 
we would need to intersect 250 non-convex regions, which is an extremely hard 
computational problem (see, e. g., [5]). 
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Figure 5. Extended Voronoi region over a tessellation 



To avoid this we do the following, at the expense of getting a slightly larger 
(hence, for our purposes worse) extended Voronoi region: instead of obtaining the 
extended Voronoi region directly as an intersection, we obtain it as a union of tiles 
of the same auxiliary tessellation T . That is to say: 

(1) We first identify a finite (but large) population of tiles in T that form 
themselves an extended Voronoi region (that is, which are guaranteed to 
contain Vorc p (p) for every p £ To). 

(2) We then process one by one the rotations and translations we are interested 
in, and at each step discard from our initial population the tiles that do 
not intersect the corresponding region of Lemma 12.41 or Lemma 12.51 

What we eventually obtain with this method is the set of tiles of T that intersect 
the region of Figure 01 See Figure [5] 

The reader may wonder why step (2) above is computationally not so hard, since 
it again involves the intersection of non-convex regions. The reason is that now we 
intersect only two such regions at a time (one tile and one strip or dihedron) and 
do not need to intersect the resulting set with anything else. Also (although less 
important), we do not really need to compute an intersection but only to check its 
cmptyness. 

3. Detailed discussion of the method 

3.1. The auxiliary tessellation. The first natural choice of an auxiliary tessel- 
lation would be a tessellation by fundamental domains of the normalizer Af(G), or 
of any subgroup of it. In this way all the tiles are equivalent: By Lemma 12.31 in 
such a tessellation we would need to compute only one extended Voronoi region. 

Since we have Af(Q) < A/"(G) for every quarter group, as a first step let us 
compute a fundamental domain of J\f(Q) (later, we will subdivide it further). In 
order to get a fundamental domain with some symmetries, we start with the Voronoi 
region of one of its degenerate orbits, namely the orbit with base point at the origin 
(0,0,0). This orbit is a body centered cubic lattice of side 1/2. Its Voronoi cell is 



DIRICHLET STEREOHEDRA FOR QUARTER CUBIC GROUPS 



9 



Figure 6. A truncated octahedron, the Voronoi cell of the body 
centered cubic lattice 



a truncated octahedron (see Figure EJ) whose 14 facets have supporting inequalities 
111 3 

±X < — , ±W < — , ±Z < — , ±X ± H ± Z < — . 

- 4' y - 4' - 4' y - 8 

The 24 vertices of the Voronoi cell are 

(±±±±0), (±io,±i), (o,±±,±§), 

(±i,±i,0), (±io,±i), (0,±±,±±). 

The stabilizer of this truncated octahedron inAf(Q) has order six. It is generated 
by the central symmetry (inversion) around the origin and the rotation of order 
three with axis in the line x = z = y. Hence, one way of obtaining a fundamental 
domain for Af(Q) is to intersect the truncated octahedron with a sector of angle 
7r/3 with edge in that line. There are several inequivalent ways of doing it, but we 
choose the sector defined by the inequalities x < y < z. 

In this way, we obtain as fundamental domain of N(Q) the polytope with the 
following 12 vertices: 

(1 I I) 

V8> 8' 8/ 

V4> 16' 16/ \4> 8' / Vl6'16' / 

(I =1=1) (I -±) (i =±) (i J- ^1 

V4' 16 ' 16 / V4' U J 8/ V8 ' u ' 4 / V 16 ' 16' 4 / 

V u ' 16 ' 16 / V u ' 8 ' 4 / V 16 ' 16 ' 4 / 



This polytope is depicted on the left side of Figure [7] It has eight facets: the two 
"big ones" in the planes x = y and y = z are pentagons and the other six are parts of 
facets of the truncated octahedron: two halves-of-squares, two halves-of-hexagons 
and two sixths-of-hexagons. 

As we said above, to get better bounds, we further subdivide this fundamental 
domain of Af(Q). We do it cutting with the coordinate planes z — 0, y — and 
x = 0. We get the following four sub-polytopes Aq, Bo, Cq and Dq (see Figure [7]): 

. o=c ^ (o , o , o) ,(-,0,(-,^),(-, o ),(3,3 iO Wl Oio)) 
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Figure 7. Two views of the fundamental domain of Af(Q) and its 
subdivision into four prototiles, Aq, Bq, Cq and Dq 



conv ( (0,0,0), [ j, 0,0 



,0, 



r:0, 



11-1 

16' 16'~4~ 



0,0, 



Co := co„v((0,0,0),(i,0,0),(i,^4),(i,04), 



conv (0, 0, 0) , 



0, 



3 -3 



-1 -1 



8 ' 4 



-1 -1 -1 



-1 -1 -1 



16 ' 16 J ' V 16 ' 16 ' 4 



0,0, 



1 



Definition 3.1. From now on, the tessellation 7~ o/R 3 obtained letting J\f(Q) act 
on the four "prototiles" Aq, Bq, Cq, and Dq is called auxiliary tessellation. 



3.2. The initial population of tiles. As said in Scction[2j4j we compute extended 
Voronoi regions starting with a finite number of tiles of T whose union in guaranteed 
to contain Vorc p (p) for every p £ To, and then eliminating as many of them as we 
can. As initial population of tiles, we take all the ones contained in the cube [—1, l] 3 . 
That this is enough follows from the fact that the four prototiles are contained in 
the cube [—1/4, 1/4] 3 and our group Q contains the translations of length 1 in the 
three coordinate directions. Lemma 12.51 then implies that the extended Voronoi 
region (of any of the prototiles) is contained in the cube [—3/4, 3/4] 3 . 

Since the four prototiles together form a fundamental domain of J\f(Q) and since 
the density of (any generic orbit of) J\f(Q) is 96 points per unit cube, in the initial 
population we have 96 x 8 = 768 tiles of each of the four types, and 3 072 tiles in 
total. 

In our implementation, each tile of T is characterized by its "type" (A, B, C or 
D) and the transformation p 6 Af(Q) that sends the corresponding prototile (Aq, 
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1/2 
1/2 


1 






-1 








1/2 


1 






1 






1 -1 
1 




-1 








1/2 








1 











1 


1/2 








-1 


1 





-1 




1 











1 











1 











1 











Cs 


3/4 





-1 


1/4 





1 





3/4 





1 





1/4 





-1 





1/4 


-1 








1/4 


1 








3/4 


-1 








3/4 


1 










3/4 








1 


1/4 








1 


1/4 








-1 


3/4 








-1 



Table 3. The 32 matrices of the isometrics which send horizontal 



planes to horizontal planes for Ao 



Bo, Co or Do) to it. We compute the initial tiles of type, say, A, in the initial 
tessellation as follows: 

• From the graphic representation of M(Q) in Figurc[2]we read the 32 isome- 
tries in Af(Q) that send horizontal planes to horizontal planes and send a 
base point p € [0, 1/4] 3 to lie in the unit cube [0, l] 3 . The matrices of these 
isometries are shown in Table [3] For future reference, these 32 matrices 
appear in eight subsets Ci, i = 1, . . . , 8, representing the cosets modulo Q. 
Table [4] shows which cosets of isometries form each of the eight quarter 
groups. 

• To these 32 isometries we apply the triad rotation with axis x = y = z, to 
obtain the 96 elements of M(Q) that send a point p e A C [0, 1/4] 3 to lie 
in [0, l] 3 . 

• Finally, we compose these 96 with the translations by the vectors (0,0,0), 
(-1,0,0), (0,-1,0), (0,0,-1), (-1,-1,0), (-1,0,-1), (0,-1,-1) and 
(—1, —1, —1) in order to get all the tiles of type A contained in [—1, l] 3 . 

With the tiles of the other types, B, C and D, we work the same way, except 
for the fact that the prototiles are now not in [0, l] 3 . For example: Bo is in the 
unit cube [0, 1] x [0, 1] x [—1, 0]. Still, to build our initial population of B-tiles, we 
want to start with the 32 horizontal transformations that send Bq to lie in [0, l] 3 . 
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Quarter groups 


Cosets 


q a 


Ci , . . . , Cg 


74i32 


Ci , C% , C3 , C4 


743d 


Ci, C3, Ce, Cg 


'P 


Ci , C3 , C5 , C7 


P4i32 


Ci, C2 


J2'3 


Cl,C3 


P^3 

a 


Cl,C7 


P2 1 3 


Ci 



Table 4. Cosets of isometries in each quarter group 





Type 5 


Ci 


10 
10 
10 
10 1 




10 
1/2 -10 
1/2 10 

-1 




c 3 




10 
10 
1/2 0-10 
-1 




10 
1/2 -10 
1 0-10 
10 1 


c 5 


10 
1/2 10 
10 
-1 




1 
1 -10 
1/2 10 
1 1 




c 7 




10 
1/2 10 
1/2 0-10 

10 1 




10 
1-10 
10-10 
-1 



Table 5. The new matrices of the isometries which send horizon- 



tal planes to horizontal planes for Bq 



These are the same 32 of Table [3] except some of them need to be composed with 
a suitable integer translation. The matrices that need modification are shown in 
Tables [SJ El and [3 for types B, C and D, respectively. As in case A, once we have 
got these isometries, we apply to them the triad rotations with axis x = y = z and 
the translations by the vectors (0,0,0), (-1,0,0), (0,-1,0), (0,0,-1), (-1,-1,0), 
( — 1, 0, —1), (0, —1, —1) and ( — 1, —1, —1) to get an initial population of 96 x 8 = 768 
tiles of each type. 

3.3. Cutting with translations and rotations. At this point we can say that 
we have already computed a single extended Voronoi region for the eight groups 
and the four prototilcs: the cube [—1, l] 3 decomposed into 768 x 4 = 3072 tiles, 
768 of each type A, B, C or D. But this region is certainly not good enough for 
our purposes. Our next step is to apply Lemmas 12.41 and 12.51 to exclude from this 
initial population as many tiles as we can. 

That is to say: Let G be one of the eight quarter groups and let To be one of the 
four prototiles Aq, Bo, Co and Do- For each choice of G and To we do the following: 
• Let S be a set of rotations and translations present in G. Below we specify 
our choice of 5" for each G, but essentially what we do is we take the 
translations of smaller length and the rotations of axis closer to the origin. 
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Type C 




1 



















1 











1 











(1, 
w 





1 
















1/2 


-1 








1/2 


1 








1 





1 













1/2 





1 











-1 







1 








1 



















-1 


1/2 








-1 




1 











1 



















1 











c 3 


1 


-1 











1 
















1/2 


-1 








1 





1 





1/2 





-1 



















-1 







1/2 








-1 











-1 










1 








1 




1 



















1 











1 











c 5 


1/2 


1 
















1 


-1 











1 








1 





1 













1/2 





1 











-1 
















-1 










1 








1 


1/2 








1 




1 











1 



















1 











c 7 


1/2 
1 


-1 






1 






1/2 
1/2 


1 






-1 














1 - 



-1 






-1 








1/2 








1 


1 








1 



















-1 



Table 6. The new matrices of the isometries which send horizon- 



tal planes to horizontal planes for Co 





Type D 




1 











1 











1 











1 













1 


1 











-1 








1/2 


-1 








1/2 


1 








1 





1 





1/2 





-1 





1/2 





1 











-1 







1 








1 


1/2 








1 











-1 


1/2 








-1 




1 











1 



















1 
















-1 








1 


1 
















1/2 


-1 








1 





1 





1/2 





-1 



















-1 







1/2 








-1 











-1 










1 








1 




1 



















1 











1 











c 5 


1/2 


1 



















-1 








1 


1 








1 





1 













1/2 





1 











-1 
















-1 










1 








1 


1/2 








1 




1 











1 











1 











1 











c 7 


1/2 
1 


-1 






1 






1/2 
1/2 


1 






-1 






1 

1/2 


1 






1 






- 



-1 






-1 








1/2 








1 


1 








1 


1/2 








-1 











-1 



Table 7. The new matrices of the isometries which send horizon- 



tal planes to horizontal planes for Dq 



• For each of the 768 x 4 tiles in the initial population and for each element 
p G S, apply Lemma 12^41 (if p is a rotation) or 12.51 Tif p is a translation) and, 
if the tile is fully contained in the region of space forbidden by the Lemma, 
then remove this tile from the extended Voronoi region. 

There is a subtle point in this latter test. Let T be a tile to which we can 
apply the test of Lemma, sav. 12.41 The "forbidden region" is the unbounded 
dihedron defined by two half-spaces H[ and H'^ (see Figure [3]). What we 
test is whether each of the vertices of T belongs to each of the two half- 
spaces. But the decision of whether the tile is discarded or not depends on 
whether the dihedron is convex or not. That is, whether the angle a of the 



14 
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FRANCISCO SANTOS 



Sets 


Isommetries 


ft 


Coordinate translations: 
(±1,0,0), (0,±1,0),(0,0,±1) 


s 2 


Diagonal translations: 
(±1/2, ±1/2, ±1/2) 


s 3 


Triad rotations (Table EJ) 


($4 


Diad rotations (Table H01 


s 5 


Diad rotations (Table 1111) 


Se 


Diad rotations (Table [121 



Group 


Transl. 


Rotations 


P2i3 


ft 


ft 


P4i32 


ft 


ft, ft 


J43d 

7 F 
J2'3 


ft, ft 


ft, ft 


g d 

74i32 


ft, ft 


ft, ft, ft 



Table 8. Transformations that we use in each quarter group 



rotation p plus the angle with which To is seen from I is bigger or smaller 
than 180 degrees: 

— If this angle is bigger than 180 degrees (that is, if the forbidden region 
is convex) we discard the tile T whenever "both half-spaces H[ and 
H' 2 contain all the vertices of T" . 

— If this angle is smaller than 180 degrees (that is, if the forbidden region 
is not convex) we discard the tile T whenever "one of the half-spaces 
H[ and H 2 contains all the vertices of T" . 

If p is a translation, the second rule is always applied. 

The particular translations and rotations that we use are the following ones. 
Tables |S| El EE EH and [TJ list them more explicitly. 

• The integer translations of length one appear in all the quarter groups. In 
the five groups with a body centered translational subgroup we also use the 
half-integer diagonal translations (±1/2, ±1/2, ±1/2). 

• The triad rotations also appear in all the quarter groups. The set S3 of 
triad rotations that we use are listed in Table [9l We have included the 
ones with axis closer to the prototile, since the ones that are not close will 
cut out only portions of the extended Voronoi region that were already cut 
by other, closer rotations. The decision on where to put the threshold is 
somehow heuristic but should not affect the result much. Actually, many 
of the rotations listed in our table turned out to be already superfluous. 

• Table[TU]lists the set S4 of diad rotations with axes parallel to the coordinate 
axes that we use in the groups J^3|, J4i32, 743d, 7|3 and 72'3. 

• Table QT| lists the set S5 of diad rotations with axes parallel to the diagonals 
of the faces of the unit cube that we use in the groups ^^"3^ and P4i32. 
In the group P4i32 only some of these diad rotations are present. The list 
S 6 of them is in Table [TU 

In these four tables (Tables El [TU1 [TT1 andll2 p . the two entries in each column are 
a point and the direction of the rotation axis. 

3.4. Planar projection. In the extended Voronoi region computed so far only 
rotations and translations present in our group are taken into account. If we wish 
to deal with other transformations (glide reflections, screw rotations, etc) we could 
state the analogues of Lemmas 12.41 and 12.51 for them, but the result would be too 
complicated, perhaps impractical. 

Instead, let us consider the "horizontal subgroup" of our group G. That is, the 
subgroup G z that fixes the third coordinate. Each of the eight parts of Figure [5] 
shows the planar subgroup of the quarter group that occupies the same position in 
Figure [2] Only four different planar subgroups arise, of types pi, p2, pg and pgg. 
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Triad rotations 



Point 


Vector 


Point 


Vector 


Point 


Vector 


(0,0,0) 


(1.1.1) 


(-1/2,0,0) 


(-1, 


-1,1) 


(0,1/2,0) 


(-1 


1 


1) 


(-1,0,0) 


(1.1.1) 


(1/2,0,0) 


(-1, 


-1,1) 


(0,-1/2,0) 


(-1 


1 


1) 


(1,0,0) 


(1,1,1) 


(1/2,1,0) 


(-1, 


-1,1) 


(-1,1/2,0) 


(-1 


1 


1) 


(0,-1,0) 


(1,1,1) 


(-1/2,1,0) 


(-1, 


-1,1) 


(-1,-1/2,0) 


(-1 


1 


1) 


(0,1,0) 


(1,1,1) 


(1/2,-1,0) 


(-1, 


-1,1) 


(1,1/2,0) 


(-1 


1 


1) 


(0,0,-1) 


(1,1,1) 


(-1/2,-1,0) 


(-1, 


-1,1) 


(1,-1/2,0) 


(-1 


1 


1) 


(0,0,1) 


(1,1,1) 


(-3/2,-1,0) 


(-1, 


-1,1) 


(-1,3/2,0) 


(-1 


1 


1) 


(0,0,1/2) 


(1,-1,1) 


(1/2,0,1) 


(-1, 


-1,1) 


(0,1/2,-1) 


(-1 


1 


1) 


(1,0,1/2) 


(1,-1,1) 


(-1/2,0,1) 


(-1, 


-1,1) 


(0,-1/2,-1) 


(-1 


1 


1) 


(0,0,-1/2) 


(1,-1,1) 


(1/2,0,-1) 


(-1, 


-1,1) 


(0,1/2,1) 


(-1 


1 


1) 


(1,0,-1/2) 


(1,-1,1) 


(-1/2,0,-1) 


(-1, 


-1,1) 


(0,-1/2,1) 


(-1 


1 


1) 


(0,1,-1/2) 


(1,-1,1) 


(0,0,3/2) 


(1,- 


-1,1) 


(1,-1/2,0) 


(-1 


1 


1) 


(0,1,1/2) 


(1,-1,1) 


(0,0,-3/2) 


(1,- 


-1,1) 


(1,1/2,0) 


(-1 


1 


1) 






(1/2,3/2,0) 


(1,- 


-1,1) 


(-1,1/2,0) 


(-1 


1 


1) 



Table 9. A set 53 of triad rotations in Q and, therefore, in all 
the quarter groups 



Diad rotations parallel to the coordinate axes 



Point 


Vector 


Point 


Vector 


Point 


Vector 


(0, -3/4, -1/2) 


(1,0,0) 


(-1/2,0,-3/4) 


(0,1,0) 


(-3/4,-1/2,0) 


(0,0,1) 


(0,-3/4,0) 


(1,0,0) 


(-1/2,0,-1/4) 


(0,1,0) 


(-3/4,0,0) 


(0,0,1) 


(0,-3/4,1/2) 


(1,0,0) 


(-1/2,0,1/4) 


(0,1,0) 


(-3/4,1/2,0) 


(0,0,1) 


(0,-1/4,-1/2) 


(1,0,0) 


(-1/2,0,3/4) 


(0,1,0) 


(-1/4,-1/2,0) 


(0,0,1) 


(0,-1/4,0) 


(1,0,0) 


(0,0,-3/4) 


(0,1,0) 


(-1/4,0,0) 


(0,0,1) 


(0,-1/4,1/2) 


(1,0,0) 


(0,0,-1/4) 


(0,1,0) 


(-1/4,1/2,0) 


(0,0,1) 


(0,1/4,-1/2) 


(1,0,0) 


(0,0,1/4) 


(0,1,0) 


(1/4,-1/2,0) 


(0,0,1) 


(0,1/4,0) 


(1,0,0) 


(0,0,3/4) 


(0,1,0) 


(1/4,0,0) 


(0,0,1) 


(0,1/4,1/2) 


(1,0,0) 


(1/2,0,-3/4) 


(0,1,0) 


(1/4,1/2,0) 


(0,0,1) 


(0,3/4,-1/2) 


(1,0,0) 


(1/2,0,-1/4) 


(0,1,0) 


(3/4,-1/2,0) 


(0,0,1) 


(0,3/4,0) 


(1,0,0) 


(1/2,0,1/4) 


(0,1,0) 


(3/4,0,0) 


(0,0,1) 


(0,3/4,1/2) 


(1,0,0) 


(1/2,0,3/4) 


(0,1,0) 


(3/4,1/2,0) 


(0,0,1) 



Table 10. A set 54 of diad rotations parallel to the coordinate 
axes that appear in the groups /— 3|, I4i32, /43d, 7-3 and J2'3 



As it turns out, extended Voronoi regions (and influence regions) for these four 
groups were computed in [3] and [T] . Our idea is to take advantage of the following 
simple fact: 

Lemma 3.2. Let ir : R 3 — > R 2 denote the vertical projection, omitting the third 
coordinate. Let Tq be a prototile and let G be a crystallographic group. Let G z 
denote the horizontal subgroup of G, that is the set of transformations that send 
each horizontal plane to itself. Then: 

7r(VorExt G (T )) C VorExt Gz (tt(T )). 

Proof. The only thing to prove is the same formula for Voronoi regions. That is, 
for every base point p 6 To we need to show that: 

7r(Vor Gp (p)) C Voi- Gz{7i{p)) (tt(p)). 
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Point 


Vector 


Point 


Vector 


Point 


(-7/4,0, 1/8) 


(-1, 1,0) 


(-5/4, 0, 3/8) 


(1,1 


0) 


(1/8,-7/4,0) 


(-7/4, 0, 5/8) 


(-1, 1,0) 


(-5/4, 0, 7/8) 


(1,1 


0) 


(5/8,-7/4,0) 


(-5/4,0,3/8) 


(-1, 1,0) 


(-3/4, 0, 1/8) 


(1,1 


0) 


(3/8,-5/4,0) 


(-5/4, 0, 7/8) 


(-1, 1,0) 


(-3/4, 0, 5/8) 


(1,1 


0) 


(7/8,-5/4,0) 


(-3/4,0,1/8) 


(-1, 1,0) 


(-1/4,0,3/8) 


(1,1 


0) 


(1/8,-3/4,0) 


(-3/4, 0, 5/8) 


(-1, 1,0) 


(-1/4,0,7/8) 


(1,1 


0) 


(5/8,-3/4,0) 


(-1/4,0,3/8) 


(-1, 1,0) 


(1/4,0,1/8) 


(1,1 


0) 


(3/8, -1/4, 0) 


(-1/4,0,7/8) 


(-1, 1,0) 


(1/4,0,5/8) 


(1,1 


0) 


(7/8,-1/4,0) 


(1/4,0,1/8) 


(-1,1,0) 


(3/4, 0, 3/8) 


(1,1 


0) 


(1/8, 1/4,0) 


(1/4,0,5/8) 


(-1, 1,0) 


(3/4, 0, 7/8) 


(1,1 


0) 


(5/8, 1/4, 0) 


(3/4,0,3/8) 


(-1, 1,0) 


(5/4,0,1/8) 


(1,1 


0) 


(3/8, 3/4, 0) 


(3/4,0,7/8) 


(-1, 1,0) 


(5/4,0,5/8) 


(1,1 


0) 


(7/8, 3/4, 0) 


(5/4,0,1/8) 


(-1, 1,0) 


(7/4, 0, 3/8) 


(1,1 


0) 


(1/8, 5/4, 0) 


(5/4,0,5/8) 


(-1, 1,0) 


(7/4, 0, 7/8) 


(1,1 


0) 


(5/8, 5/4, 0) 


(7/4,0,3/8) 


(-1, 1,0) 


(-5/4,0,-1/8) 


(1,1 


0) 


(3/8, 7/4, 0) 


(7/4,0,7/8) 


(-1, 1,0) 


(-5/4,0,-5/8) 


(1,1 


0) 


(7/8, 7/4, 0) 


(-7/4,0,-3/8) 


(-1, 1,0) 


(-3/4,0, -3/8) 


(1,1 


0) 


(-3/8, -7/4,0) 


(-7/4, 0, -7/8) 


(-1, 1,0) 


(-3/4, 0, -7/8) 


(1,1 


0) 


(-7/8, -7/4,0) 


(-5/4,0,-1/8) 


(-1, 1,0) 


(-1/4,0,-1/8) 


(1,1 


0) 


(-1/8, -5/4,0) 


(-5/4,0,-5/8) 


(-1, 1,0) 


(-1/4,0,-5/8) 


(1,1 


0) 


(-5/8, -5/4,0) 


(-3/4, 0, -3/8) 


(-1, 1,0) 


(1/4,0,-3/8) 


(1,1 


0) 


(-3/8, -3/4,0) 


(-3/4, 0, -7/8) 


(-1, 1,0) 


(1/4,0,-7/8) 


(1,1 


0) 


(-7/8, -3/4,0) 


(-1/4,0,-1/8) 


(-1, 1,0) 


(3/4,0,-1/8) 


(1,1 


0) 


(-1/8, -1/4,0) 


(-1/4,0,-5/8) 


(-1, 1,0) 


(3/4, 0, -5/8) 


(1,1 


0) 


(-5/8, -1/4,0) 


(1/4,0, -3/8) 


(-1, 1,0) 


(5/4, 0, -3/8) 


(1,1 


0) 


(-3/8, 1/4, 0) 


(1/4,0,-7/8) 


(-1, 1,0) 


(5/4, 0, -7/8) 


(1,1 


0) 


(-7/8, 1/4, 0) 


(3/4,0,-1/8) 


(-1, 1,0) 


(7/4,0,-1/8) 


(1,1 


0) 


(-1/8,3/4,0) 


(3/4,0,-5/8) 


(-1, 1,0) 


(7/4, 0, -5/8) 


(1,1 


0) 


(-5/8, 3/4, 0) 


(5/4,0,-3/8) 


(-1, 1,0) 


(0, 1/8, -7/4) 


(1,0, 


-1) 


(-3/8, 5/4, 0) 


(5/4,0,-7/8) 


(-1, 1,0) 


(0, 1/8, -7/4) 


(1,0, 


-1) 


(-7/8, 5/4, 0) 


(7/4,0,-1/8) 


(-1, 1,0) 


(0,3/8,-5/4) 


(1,0, 


-1) 


(-1/8, 7/4,0) 


(7/4,0,-5/8) 


(-1, 1,0) 


(0, 7/8, -5/4) 


(1,0, 


-1) 


(-5/8, 7/4, 0) 


(3/8,-5/4,0) 


(0,1,1) 


(0, 1/8, -3/4) 


(1,0, 


-1) 


(0,3/8, -5/4) 


(7/8, -5/4,0) 


(0,1,1) 


(0, 5/8, -3/4) 


(1,0, 


-1) 


(0, 7/8, -5/4) 


(1/8, -3/4,0) 


(0,1,1) 


(0,3/8,-1/4) 


(1,0, 


-1) 


(0, 1/8, -3/4) 


(5/8, -3/4, 0) 


(0,1,1) 


(0,7/8,-1/4) 


(1,0, 


-1) 


(0, 5/8, -3/4) 


(3/8,-1/4,0) 


(0,1,1) 


(0, 1/8, 1/4) 


(1,0, 


-1) 


(0,3/8, -1/4) 


(7/8,-1/4,0) 


(0,1,1) 


(0,5/8,1/4) 


(1,0, 


-1) 


(0,7/8,-1/4) 


(1/8, 1/4,0) 


(0,1,1) 


(0,3/8,3/4) 


(1,0, 


-1) 


(0, 1/8, 1/4) 


(5/8,1/4,0) 


(0,1,1) 


(0, 7/8, 3/4) 


(1,0, 


-1) 


(0, 5/8, 1/4) 


(3/8,3/4,0) 


(0,1,1) 


(0,1/8,5/4) 


(1,0, 


-1) 


(0, 3/8, 3/4) 


(7/8,3/4,0) 


(0,1,1) 


(0,5/8,5/4) 


(1,0, 


-1) 


(0, 7/8, 3/4) 


(1/8,5/4,0) 


(0,1,1) 


(0,3/8,7/4) 


(1,0, 


-1) 


(0, 1/8, 5/4) 


(5/8,5/4,0) 


(0,1,1) 


(0, 7/8, 7/4) 


(1,0, 


-1) 


(0, 5/8, 5/4) 


(3/8,7/4,0) 


(0,1,1) 


(0, -3/8, -7/4) 


(1,0, 


-1) 


(0, 3/8, 7/4) 


(7/8,7/4,0) 


(0,1,1) 


(0, -7/8, -7/4) 


(1,0, 


-1) 


(0, 7/8, 7/4) 


(-1/8,-5/4,0) 


(0,1,1) 


(0,-1/8, -5/4,0) 


(1,0, 


-1) 


(0, -1/8, -5/4) 


(-5/8, -5/4,0) 


(0,1,1) 


(0, -5/8, -5/4) 


(1,0, 


-1) 


(0, -5/8, -5/4) 


(-3/8, -3/4,0) 


(0,1,1) 


(0, -3/8, -3/4) 


(1,0, 


-1) 


(0, -3/8, -3/4) 


(-7/8, -3/4,0) 


(0,1,1) 


(0, -7/8, -3/4) 


(1,0, 


-1) 


(0, -7/8, -3/4) 


(-1/8,-1/4,0) 


(0,1,1) 


(0, -1/8, -1/4) 


(1,0, 


-1) 


(0, -1/8, -1/4) 


(-5/8,-1/4,0) 


(0,1,1) 


(0, -5/8, -1/4) 


(1,0, 


-1) 


(0, -5/8, -1/4) 


(-3/8,1/4,0) 


(0,1,1) 


(0, -3/8, 1/4) 


(1,0, 


-1) 


(0, -3/8, 1/4) 


(-7/8,1/4,0) 


(0,1,1) 


(0, -7/8, 1/4) 


(1,0, 


-1) 


(0, -7/8, 1/4) 


(-1/8,3/4,0) 


(0,1,1) 


(0,-1/8,3/4) 


(1,0, 


-1) 


(0,-1/8,3/4) 


(-5/8,3/4,0) 


(0,1,1) 


(0, -5/8,3/4) 


(1,0, 


-1) 


(0, -5/8,3/4) 


(-3/8,5/4,0) 


(0,1,1) 


(0,-3/8,5/4) 


(1,0, 


-1) 


(0,-3/8,5/4) 


(-7/8,5/4,0) 


(0,1,1) 


(0, -7/8, 5/4) 


(1,0, 


-1) 


(0, -7/8, 5/4) 


(-1/8,7/4,0) 


(0,1,1) 


(0,-1/8,7/4) 


(1,0, 


-1) 


(0,-1/8,7/4) 


(-5/8,7/4,0) 


(0,1,1) 


(0,-5/8,7/4) 


(1,0, 


-1) 


(0,-5/8,7/4) 



Vector 



(0, 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0, 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0, 
(0. 
(0. 
(0. 
(0. 



(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 


(1 


0. 


1) 



Table 11. A set S$ of diad rotations parallel to the diagonal of 
the faces of the cube, in the groups and I4i32 
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Diad rotations parallel to the diagonals of the faces of the unit cube 



Point 



Vector 



Point 



Vector 



Point 



Vector 



(-5/4,0,3/8) 
(-5/4, 0, 7/8) 
(-1/4,0,3/8) 
(-1/4,0,7/8) 
(3/4,0,3/8) 
(3/4,0,7/8) 
(7/4,0,3/8) 
(7/4,0,7/8) 
(-5/4,0,-1/8) 
(-5/4,0, -5/8) 
(-1/4,0,-1/8) 
(-1/4,0,-5/8) 
(3/4,0,-1/8) 
(3/4,0, -5/8) 
(7/4,0,-1/8) 
(7/4,0,-5/8) 
(1/8,-3/4,0) 
(5/8,-3/4,0) 
(1/8,1/4,0) 
(5/8,1/4,0) 
(1/8,5/4,0) 
(5/8,5/4, 0) 
(-3/8, -3/4,0) 
(-7/8,-3/4,0) 
(-3/8,1/4,0) 
(-7/8,1/4,0) 
(-3/8, 5/4, 0) 
(-7/8,5/4,0) 



-1,1, 
-1,1, 
-1,1, 
-1,1, 
-1, 1. 

-1, 1, 
-1, 1, 
-1,1,0) 
-1,1,0) 
-1, 1, 
-1, 1, 
-1,1, 
-1,1, 
-1, 1, 
-1, 1, 
-1,1,0) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 
(0,1,1) 



.0) 
,0) 
.0) 
.0) 
.0) 
.0) 
.0) 



.0) 
.0) 
.0) 
.0) 
.0) 
.0) 



(-3/4,0, 1/8) 
(-3/4, 0, 5/8) 
(1/4,0,1/8) 
(1/4,0,5/8) 
(5/4,0,1/8) 
(5/4,0,5/8) 
(-3/4, 0, -3/8) 
(-3/4, 0, -7/8) 
(1/4,0,-3/8) 
(1/4,0,-7/8) 
(5/4, 0, -3/8) 
(5/4,0,-7/8) 
(0,3/8,-5/4) 
(0, 7/8, -5/4) 
(0,3/8,-1/4) 
(0,7/8,-1/4) 
(0,3/8,3/4) 
(0, 7/8, 3/4) 
(0, 3/8, 7/4) 
(0, 7/8, 7/4) 
(0,-1/8,-5/4,0) 
(0, -5/8, -5/4) 
(0,-1/8,3/4) 
(0, -5/8,3/4) 
(0,-1/8,-1/4) 
(0, -5/8, -1/4) 
(0,-1/8,7/4) 
(0, -5/8, 7/4) 



(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 
(1,1,0) 



(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 
(1,0, 



-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 



(3/8, -5/4, 0) 
(7/8,-5/4,0) 
(3/8,-1/4,0) 
(7/8,-1/4,0) 
(3/8,3/4,0) 
(7/8,3/4, 0) 
(3/8, 7/4, 0) 
(7/8, 7/4, 0) 
(-1/8,-5/4,0) 
(-5/8, -5/4,0) 
(-1/8,-1/4,0) 
(-5/8,-1/4,0) 
(-1/8,3/4,0) 
(-5/8, 3/4, 0) 
(-1/8, 7/4,0) 
(-5/8, 7/4, 0) 
(0, 1/8, -3/4) 
(0, 5/8, -3/4) 
(0, 1/8, 1/4) 
(0, 5/8, 1/4) 
(0, 1/8, 5/4) 
(0, 5/8, 5/4) 
(0, -3/8, -3/4) 
(0, -7/8, -3/4) 
(0, -3/8, 1/4) 
(0, -7/8, 1/4) 
(0,-3/8,5/4) 
(0,-7/8,5/4) 



(0. 
(0. 
(0, 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0. 
(0, 
(0. 



(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 
(1,0,1) 



Table 12. A set Sq C S5 of diad rotations parallel to the diagonal 
of the faces of the unit cube that appear in the group P4i32 



The latter follows from the fact that G z (tt(p)) equals the projection of the part of 
Gp in the same horizontal plane as p (for a sufficiently generic p). □ 



That is to say: we can intersect our previously computed extended Voronoi 
regions VorExt G (T ) (for T G {A , . . . , A)}) with tt' 1 (VorExt Gz (tt(T ))) and use 
the new, smaller region to compute influence regions. In order to do this we check, 
for each tile in VorExto(To), whether its projection intersects VorExtg^ (vr(To)). If 
it does not then we discard that tile. Needless to say that we apply this process 
three times to each prototile, for the three coordinate projections XY, XZ and 
YZ. 

For the groups with a planar pi this step does not reduce the regions at all, so 
we do not show the details. But for the other three types of groups, p2, pg and pgg, 
this step is significant. Roughly speaking, it allows us to use some glide reflections 
present in some of the quarter groups to cut the extended Voronoi region further. 

As a planar prototile we have used a square of side 1/4 for pg and half of it 
(a right triangle) for pgg and p2. The reason for this difference is that in pg the 
extended Voronoi region would not decrease significantly if we took a triangle. We 
omit the calculation of the extended Voronoi regions, shown in Figure [5] The 
details can be found in [3] and [1] . Let us only say that the calculation is based on 
the planar version of Lemmas 12.41 and 12.51 plus the analogue for glide reflections. 
Roughly speaking, each glide reflection with translation vector, say v, "acts as if" 
the group G z contained a rotation of order two with center on the glide reflection 
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Figure 8. Plane groups of the quarter groups. 



axis and at distance \v/2\ from the orthogonal projection of the prototile To to the 
axis. Glide reflections are indicated as dotted lines in the figures. 

Of course, we need to check whether each prototile T that we used in the 3D 
step projects to lie in one planar tile. That is, in one tile obtained from the one 
in Figure by the action of the planar group G z . If this is not the case then the 
planar extended Voronoi region we consider is going to be a union of more than 
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(a) plane group pg. (b) plane group pgg. 



(c) plane group p2. 



Figure 9. Extended Voronoi Regions of the plane groups pg, pgg 
and p2. 



one of the regions shown in the figure. The result of this check, for each of the four 
tiles Aq, Bq, Cq and Dq, each of the three coordinate projections, and the three 
possible plane groups, is illustrated in Tables ITOl and ITU As shown in the pictures, 
for the group pg one planar tile is always enough, but for the groups pgg and p2 
we usually need two. 

3.5. Influence regions. By Theorem 12. 21 once we know the influence region of, 
say, the prototile Aq, for a group G, we can obtain an upper bound for the maximum 
number of facets of the Dirichlct stereohedra of G with base point in Aq. We only 
need to count how many orbit points lie in the influence region, that is, how many 
tiles in the G-orbit of Aq are listed in our influence region. (Remember that, by 
construction, our influence region is a union of tiles of the auxiliary tessellation T.) 
Let us see in more detail how we compute this number. 

By definition, a tile T is in Infl G (A ) if and only if VorExt G (A ) n VorExt G (T) ^ 
0. Note that T can a priori be a tile of any of the four types, A, B, C or D. But 
for T produce an orbit point it is necessary (but not sufficient, see below) that T 
be of the same type, A, as our prototile. Such a tile A' is in Infl G (Ao) if and only 
if there is a tile T" (of any of the four types) in VorExt G (A)) H VorExt G (yl'). 

At this point we have to clarify one feature of our encoding of tiles that is useful 
here. Remember that for us VorExt G (^4o) is a list of tiles of T, and that each tile 
T of T has been encoded as the transformation p € Af(Q) that sends the prototile 
To of the same type to T (see Figure [TO]). Suppose now that A' is a tile of type A 
in Infl G (Ao). That is, A is such that 

T' £ VorExt G (A>) n VorExt G (A') 

for some T'. Let fi G Af(Q) be the transformation that sends Aq to A'. Then, by 
Lemma [23] (and since Af(Q) C Af(G) for every quarter group G), VorExt G (A') = 
p, VorExt G (A). That is, there is a second tile T\ in VorExt G (A) such that T" = fxTi. 
Let p and p' denote the transformations that send the prototile To of the appropriate 
type to T' and T\ , respectively. We obviously have that (see Figure [TO] again) : 

M = P °P~ l - 

Therefore, all the tiles A C Infl G (^4o) that may possibly produce a facet in our 
Dirichlet stereohedra can be represented as the composition p 1 o where p and 
p' are two isometries (that represent tiles) in VorExt G (^4o)- 

There is a final step. Each transformation p = p' o p^ 1 obtained in this way is 
clearly in Af(Q), but it may not be in our particular group G. Only if it is in G it 
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Plane XY 


Plane XZ 


Plane YZ 



Table 13. Extended Voronoi Regions of the plane groups of the 
quarter groups. Prototiles of types A and B 



contributes one neighbor to the bound. To check this we use the coset classification 
of all the tiles in our initial population, shown in Table 0J 
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Plane XY 



Plane XZ 



Plane YZ 



Table 14. Extended Voronoi Regions of the plane groups of the 
quarter groups. Prototiles of types C and D 
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Figure 10. Calculation of the influence region 



The influence regions for the three other prototiles Bq, Co and Dq are obtained 
in exactly the same way. Once we have the number of elements in each influence 
region of G for each prototile, we take the biggest of these numbers as a bound of 
the maximum number of facets of Dirichlet stereohedra for G. 

The results of these calculations are shown in columns (1), (2), (3) and (4) 
of Table [T] For each group, the rightmost number of these three columns is the 
maximum size of influence regions obtained. Column (1) shows the bound obtained 
if we neglect the diad rotations and the planar projection step. Columns (2), (3) and 
(4) show how the bound decreases after considering, respectively, the coordinate 
diad rotations, the diagonal diad rotations, and the planar projection step. 
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